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Ballistic transport through a collection of quantum billiards in undoped graphene is studied 
analytically within the conformal mapping technique. The billiards show pseudodiffusive behavior, 
with the conductance equal to that of a classical conductor characterized by the conductivity oq — 
4e^/7r/i, and the Fano factor F — 1/3. By shrinking at least one of the billiard openings, we observe 
a tunneling behavior, where the conductance shows a power-law decay with the system size, and 
the shot noise is Poissonian {F = 1). In the crossover region between tunneling and pseudodiffusive 
regimes, the conductance G [l — F) x se^ /h. The degeneracy s — 8 for the Corbino disk, which 
preserves the full symmetry of the Dirac equation, s = 4 for billiards bounded with smooth edges 
which break the symplectic symmetry, and s — 2 when abrupt edges lead to strong intervalley 
scattering. An alternative, analytical or numerical technique, is utilized for each of the billiards to 
confirm the applicability of the conformal mapping for various boundary conditions. 

PACS numbers: 73.50.Td, 73. 23. Ad, 73.63.-b 



I. INTRODUCTION 

The isolation of single layers of carbon (graphene) 
whose low-energy spectrum is described by the Dirac- 
Weyl Hamiltonian of massless spin- 1/2 fermions [V, has 
offered physicists the unique possibility to test the pre- 
dictions of relativistic quantum mechanics in a condensed 
phase. A particular attention focuses on ballistic trans- 
port [2, , as the unusual band structure of a carbon mono- 
layer 1^ leads simultaneously to a divergent Fermi wave- 
length Xp ^ oo in the undoped graphene limit, and 
to a zero bandgap. For these reasons, the quantum- 
mechanical wave character of an electron plays an es- 
sential role in transport even through a macroscopic 
graphene sample, provided that the influence of disor- 
der is negligible A separate issue concerns the fact 
that Dirac fermions in graphene occur in two degener- 
ate families, resulting from the presence of two differ- 
ent valleys in the band structure. This valley degree of 
freedom offers conceptually new possibilities to control 
charge carriers — the so-called "valleytronics" 0. 

So far, extensive theoretical studies of ballistic trans- 
port, based on mode-matching analysis for the Dirac 
equation [BJ [7] , are available for a rectangular graphene 
sample of width W, length C, and various types of 
boundary conditions. In the regime of large aspect ra- 
tios yV/£ ^ 1, the conductance of an undoped sam- 
ple scales as G = aoW/C, with the universal conductiv- 
ity (To = ie'^/nh, regardless of boundary conditions [H]. 
Moreover, as shown by Tworzydlo et al. [7], the Fano 
factor in this case coincides with that of a diffusive wire 
{F ~ 1/3). Also, the transmission eigenvalues of these 
two systems display the same distribution. This anal- 
ogy coined the term of pseudodiffusive transport, which 
describes ballistic graphene properties in the universal 
conductivity limit. 



Recent experiments report an agreement with the the- 
oretical predictions of Refs. |6J |7j for either the conduc- 
tance [S] or the Fano factor [TU]. Furthermore, the tem- 
perature dependence of the conductivity !llj also shows 
an approximate agreement with the ballistic theory gen- 
eralized to finite temperatures |12j . However, even for 
low temperatures, the convergence with W/C —^ oo is 
much slower than predicted. In particular, for the largest 
aspect ratio >V/£ = 24 studied in Ref. [TD] the deviations 
from the limiting values GC/W = ctq and F = 1/3 are 
close to 10%, whereas results of Ref. [7| show the con- 
vergence should be already reached for moderate aspect 
ratios yV/£ > 4. A clear explanation of this discrepancy 
is missing, but it is usually attributed to the fact that 
boundary conditions used in theoretical works, which de- 
scribe either an abrupt termination of a perfect honey- 
comb lattice or an infinite mass confinement |13j . may 
not model the real-sample edges correctly |14j . 

In this work, we consider ballistic graphene systems 
of geometries for which the boundary effects are absent 
or suppressed. The paper is organized as follows: In 
Sec. [n] we briefly recall the mode-matching analysis for a 
graphene strip, and show how to employ the conformal 
symmetry of the Dirac equation for undoped graphene 
|15j to obtain analytically the transmission eigenvalues 
for other systems. Then, in Sec. Ill the method is ap- 



plied to the Corbino disk. The results are compared with 
those obtained by direct mode-matching for angular mo- 
mentum eigenstates, a relation with the nonrelativistic 
electron gas in the disk setup is also discussed. In Sec. 
|IV| we study two basic billiards bounded with mass con- 
flnemcnt: a finite section of the Corbino disk and a quan- 
tum dot with circular edges. The results obtained with 
the conformal mapping technique are confirmed by the 
computer simulation of transport using the tight-binding 
model on a honeycomb lattice. Finally, in Sec. |V] we 
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study numerically the transport across an infinitely long 
nanoribbon by utilizing the 4-terminal recursive Green's 
function algorithm [TB], as well as across a finite sec- 
tion of a nanoribbon with an abrupt lattice termination. 
All the systems show pseudodiffusive transport behav- 
ior in a wide range of geometrical parameters. A fur- 
ther analogy between them appears when (at least) one 
of the leads is narrow in comparison to the character- 
istic length of the conducting region L. Namely, the 
conductance in such a quantum tunneling regime shows 
a power-law decay G oc L^"', where a is a nonuniver- 
sal (geometry-dependent) exponent. Moreover, it is re- 
lated to the shot noise by F « 1 — Gh/{se^), so that 
the Poissonian limit {F — 1) is approached for large L. 
The symmetry-dependent factor s = 8 in the presence of 
full spin, valley and symplectic degeneracy (the case of 
the Corbino disk), or s = 4 when the mass confinement 
breaks symplectic symmetry of the Dirac equation. A fi- 
nal reduction to s = 2 may be reached by adding abrupt 
(i.e. armchair) edges, which scatter the valleys. 

The original feature of the geometries studied in this 
article is that the influence of sample edges are elimi- 
nated (for the Corbino setup) or irrelevant, as the spatial 
current distribution is not uniform, but concentrated far 
away from the edges. This is why we believe that our the- 
oretical findings could be confirmed experimentally with 
better precision than that for rectangular samples, as the 
role of boundary conditions is strongly suppressed. 

II. TRANSPORT OF DIRAC FERMIONS AND 
CONFORMAL MAPPING 




FIG. 1: Quantum billiards in undoped graphene studied ana- 
lytically (schematic), (a) The Corbino disk with inner radius 
-Ri, and outer radius 7?2. (b) Generic section of the disk char- 
acterized by the spanning angle tD. (c) Quantum dot with cir- 
cular edges. A voltage source, shown on panel (a) only, drives 
the current through each of the devices. Shadow areas on all 
panels mark heavily-doped (so highly-conducting) graphene 
leads, white dots are the poles of conformal transformation 
mapping a given system onto a strip of Fig. [2] Thick lines on 
panels (b,c) indicate infinite-mass confinement. 



The compact derivation of transmission eigenvalues of 
a weakly-doped (or undoped) graphene sample coupled 
to heavily-doped graphene leads is known due to Sonin 
[17], who pointed out that one can first calculate the re- 
flection and transmission amplitudes for an interface be- 
tween weakly-doped and heavily-undoped regions, and 
then employ the double-contact formula [18]. Here we 
show that the derivation of Ref. [17j can be easily adopted 
to the Corbino disk, a finite section of it, and to a quan- 
tum dot with circular edges (all shown in Fig.jl]), as these 
systems can be obtained from a strip by appropriate con- 
formal transformations. 



A. Mode-matching for a graphene strip 



(1) (2) 




FIG. 2; Scattering on interfaces (1) and (2) between weakly- 
doped (white area) and heavily-doped (shaded area) regions 
in graphene. Horizontal dashed-dot lines mark symbolically 
generic boundary conditions applied to a strip. 



Let us first consider an electron crossing from the 
weakly-doped region {x > 0) to the heavily-doped one 
{x < 0), as depicted in Fig. [2] The Dirac Hamiltonian 
for graphene has the well-known form |19j 

Ho = vpCT ■ p, (1) 

where Wi? is the Fermi velocity, cr — {ax, cry) is the vec- 
tor operator build of Pauli matrices for the sublattice- 
pseudospin degree of freedom, and p = —ih(dx,dy) is 



the in-plane momentum operator. Due to translational 
invariance along the y-axis, the solution of the Dirac 
equation with energy E = hvpk may be written as 
^(x,?/) — xeixje'^'^yy [20], with the transverse momen- 
tum ky = Ky {ki and Ki with i = x,y denote momentum 
components in the weakly- and heavily-doped regions, re- 
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spectively), and the spinor 



-10 



-1 



1 



(2) 



iky)/k, and the hniit of an infinite 
doping (fc <C K) is imposed. The continuity of the two 
spinor components on both sides of the interface leads 
to expressions for the reflection and transmission amph- 
tudes 



- 1 



ri 



2\/cos e 



(3) 



The amphtudes ri and ti depend solely on the angle of 
incidence 9 (see Fig. [2| , what illustrates the generic fea- 
ture of ballistic transport in graphene that is insensitive 
to the lead details [H]. The reflection and transmission 
amplitudes for an electron crossing from the undoped re- 
gion to the second heavily-doped lead are r2 — r\ and 
^2 = ^1 (up to a phase factor) , as the angle of incidence 
9 —9 in this case. Thus, employing the double-contact 
formula of Ref. [IS] the total transmission probability for 
phase-coherent transport through the system of Fig. [2] is 



T - 



\tlt2\ 



|1 — 7'-^7.2e2»'/'i2 1^ 1 + (tan0sin</)i2)^ 



(4) 



where 012 = k^dx is the phase-shift earned by an elec- 
tron when passing from the first interface to the second 
one. The above expression holds true for either propa- 
gating modes (for which fc^ 



' fc2 — fc^) or, as an ana- 
lytic continuation, for evanescent modes {kx = iqx, with 
qx = \J ~ ^^)- ^'^'^ ^ confined geometry, the quantiza- 
tion of transverse momenta is determined by boundary 
conditions [6j [7] . Namely, ky — ky (with Z-integer) may 
be written in a compact form 



, ; _ gT^jl + 7) 



(5) 



where 5 = 1,2 for the closed and (generalized) periodic 
b.c, respectively. 7 = | for either mass confinement or 
antiperiodic b.c. studied in this paper. (For a nanotube- 
like geometry as considered in Ref. ^,7 = corresponds 
to periodic b.c.) 



B. Transmission via evanescent modes 

Here we limit ourselves to the case of undoped 
graphene (A; — > 0), in which the charge transport is car- 
ried fully by evanescent modes. An analytic continuation 
yields tan^ i and (f>i2 ignjC/W in Eq. Q, where 
we use the quantization (|5| and define the half-integer 
j = I + ^. As pointed out by Katsnelson and Guinea 



[15] , the zero-energy solution of the Dirac equation may 
be obtained via conformal transformation that links the 
considered geometry to a simple one, for which the wave- 
function is known |22j . In particular, if the conformal 
transformation z{w) turns the system under considera- 
tion into a rectangle of width W and length C (Fig. |2|, 
the transmission probability for the j-th evanescent mode 
may be written as 



cosh [gj In A{z(ti;)}] (Am + A-aiy 



(6) 



where j — ±i, ±7 



(with the degeneracy Tj — T-j). 



Notice that the amplitudes ([s]) remain unchanged after 
applying an arbitrary conformal transformation to the 
coordinate system of Fig. [2] so the only term in Eq. 
([4]) affected by the transformation z{w) is the phase-shift 
't'12 jfljluA. The real functional A{z{w)} is defined 

by 



lnA{z(w)} = ttC/W. 



(7) 



The explicit form of A{z{w)} depends on the geometry, 
and is given below for the examples of conformal trans- 
formation z(w) having one and two poles in a complex 
plane, which allows us to obtain expressions for transmis- 
sion probabilities through a finite section of the Corbino 
disk and through a quantum dot with circular edges, re- 
spectively. 

But first, we discuss the two basic physical regimes of 
quantum transport in graphene, which are described by 
opposite limits of Eq. ^ . The conductance of undoped 
graphene [SJ [7] is given by the Landauer formula 
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se 



E 



Tj = sTTfTo (A^^' + A- 



(8) 



with the degeneracy s = 4 (spin and valley) for smooth 
mass confinement, and s — 8 for antiperiodic b.c. due to 
an additional (symplectic) symmetry |23|. The universal 
conductivity is ctq = ie^/irh. The Fano factor also follows 
from summing over the modes 



F ■ 



_1 3 
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(9) 



J 2 ' 2 ' 



but is affected by the symmetry-dependent factors (g, s) 
only via Tj-s 

For the limit InA ^ 1, we can replace summation in 
Eq. ([8| by integration, and get 



G « Gdiff 



TTfTo 



lnA{z(«;)}' 



(10) 



where we use the relation s = Ag, valid for the two classes 
of b.c. studied here. In the In A ^ 1 limit, the relevant 



4 



information about transmission probabilities is given by 
their statistical distribution 



P{T) 



G 



diff 



(11) 



As the distribution p(T) coincides with the known dis- 
tribution p¥j for diffusion modes in a disordered metal, 
InA ^ 1 constitutes a pseudodiffusive regime of trans- 
port through graphene billiards. Notice that the generic 
conformal transformation z^w) affects p(T) only via the 
prefactor Gdiff- In particular, the Fano factor 



F= 1 - 



(T> 



1 

3' 



(12) 



regardless of the particular form of z{w). This obser- 
vation may also help to understand why experimental 
results [TD] generally show better agreement with the- 
ory for the Fano factor rather than for the conductance. 
For instance, various geometrical defects (such as a cor- 
rugation of the lead-graphene interface) may affect Gdiff 
strongly, but not affect F at all. 

In the opposite limit (In A 3> 1 ) , we find from ^ that 
Ti/2 > T3/2 > • ■ • , leading to 



G « snaoA-s, « 1 - G- 



(13) 



These expressions constitute a quantum-tunneling regime 
for ballistic graphene, in which the transport is governed 
by a single electronic mode with the four-fold (spin and 
valley) degeneracy. Below, we provide examples illus- 
trating how the power-law dependence of G on A may be 
followed by a power-law decay of G with the characteris- 
tic length-scale of the system. 




FIG. 3: Conductance and Fano factor for the undoped 
Corbino disk in graphene, as a function of the radii ratio (a,b) 
and the shot noise vs conductance diagram (c). The curves 
calculated from Eqs. Q and Q are plotted with solid lines on 
all panels. Dashed lines show the pseudodiffusive limit |17l 
on panels (a,b) and the tunneling limit (18 1 on panel (c). 



(Hereinafter, we use the symbol Log to denote the natural 
logarithm in a complex domain.) For the complex vari- 
able z — x+iy, with ^ x ^ £, and ^ y ^ W, transfor- 
mation ( [l4| leads to i?i ^ \w\ ^ R2 and ^ argw ^ 2tt 
provided the condition R2/R1 = e^'^'^l^ is satisfied. Us- 
ing ([7]), such a condition implies the functional A{z(it;)} 
to have the form 



A = A(i?i,i?2) = 



Ra 
Ri 



1/2 



(15) 



III. APPLICATION TO THE CORBINO DISK 

The Corbino setup, in which the graphene sample 
formed as an annulus is attached to coaxial leads, as 
shown schematically in Fig. [ija) , seems to be the simplest 
way to eliminate boundary effects, which are claimed to 
strongly affect experimental results for rectangular sam- 
ples with small and moderate aspect ratios [SJ [TU] . In this 
section, we first utilize the conformal mapping technique 
to find transmission eigenvalues for an undoped disk, and 
then compare the results with that obtained by a direct 
wave-function matching, possible also for a doped disk. 



As the conformal mapping is known, the only part to 
be explained now are the boundary conditions applied 
to a strip of Fig. [2j To define them, one needs to no- 
tice that after a rotation by 2tt in the coordinate sys- 
tem of Fig. [TJa), the spinor part of the wavefunction 
acquires the Berry phase [261 123 e*'^°'^ = —1. Within 
the mapping ( |14[ ), a rotation by 27r turns into a shift 
along the y-axis in Fig. [2] by a strip width W. This 
is why the spinor-rotational invariance of the original 
wavefunction implies antiperiodic boundary conditions 
^'(a;, y + W) = — ^'(a;, y) for a strip. Such boundary con- 
ditions, together with the functional A{z{w)} given by 
(151, lead the formula ^ for transmission probabilities 



to a form 



A. Conformal mapping for an undoped disk 

The conformal transformation that changes the 
Corbino disk with the inner radius i?i and the outer ra- 
dius i?2, shown in Fig.jlja), into a rectangle of the width 
yy and the length C (see Fig. |2]), is given by [25] 



w 



(14) 



cosh^ [jTn(i?i/i?2 



(16) 



A generalization for the setup with circular, but not coax- 
ial contacts is presented in Appendix |X 

The dependence of the conductance (8l and the Fano 
factor ([9]) on the radii ratio R1/R2 is plotted in Fig. [s] 
(solid lines). The limiting behavior for R1/R2 w 1, cor- 
responding In A ^ 1 (15) is characterized by G ~ Gdiff 
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(10 1, with 



G, 



diff 



27r(To 



ln(i?2/i?i)' 



and 



F 



(17) 



The formula for G^iff coincides with the well-known clas- 
sical conductance of the Corbino disk [28 . The asymp- 



totic values (17 1 are depicted with dashed red lines on 
Fig. |3];a,b). In the opposite limit (_Ri <C R2), Eq. ([Ts]) 



takes the form 



G « Sttcto 



8e2 



(18) 



The second formula from above is shown in Fig.jsjjc) with 
dashed black line. 

The results presented in Fig. [3] show that the pseu- 
dodifFusive formulas (17) for G and F match the exact 
expressions ([8| and fc| with Tj given by (16 1 in a rel- 
atively wide range ofratios i?i/i?2- Namely, the agree- 
ment becomes better than 1% if Ri ^ 0.29i?2 for the 
conductance, and if i?i ^ 0.43i?2 for the Fano factor. 
For smaller R1/R2, one can identify the crossover from 
the pseudodiffusive to quantum-tunneling behavior. In 
particular, the exact values of G are closer to the tun- 
neling formula (18 1 than to the pseudodiffusive formula 
(17 1 below R1/R2 — 0.11. The same is observed for F 
below R1/R2 = 0.16. The most characteristic feature of 
the tunneling regime is the relation G w (1 — F) x Se^/h, 



following from (18). It is satisfied with an accuracy bet- 
ter than 10% for G < 4e^/ft- (or F > 0.5), corresponding 
to R1/R2 ^ 0.2. In this range, we also find that the 
conductance decays (at fixed i?i) as G oc 1/L, where 
L = R2 — Ri « i?2 is the characteristic length of the 
sample area. 



(a) 




the mapping (Al I leads to A « Z/r and, subsequently, to 



the quadratic decay of the conductance 



G w 87rcri 



for r <ti I 



(19) 



(with the radius of each lead r and the distance between 
leads I). In the second case, the functional A w (2/'/r)^/^ 
and the conductance 



G « inao 



(20) 



shows reciprocal decay with the sample area length, sim- 
ilarly as observed for the Corbino disk. The approximate 
relation between the conductance and the Fano factor 
G w (1 — F) X Se^/h holds true for both situations of Fig. 
[4] showing the tunneling-transport regime in graphene 
appears generically for a setup consisting of (at least) 
one narrow circular lead. 



B. Electron transport at finite doping 

We complement the study of the Corbino disk in 
graphene with its transport properties at finite doping, 
characterized by the chemical potential /ig = ±hvpk 
(where fiQ > and fiQ < refers to electron and hole 
doping, respectively). The analysis is closely related to 
that for the electronic levels of graphene rings p9] . The 
single valley Hamiltonian for the doped disk reads 



H = H, 







U{r)ao, 



(21) 



where Hq is given by ([T]), the electrostatic potential 
r € i?9. and U{r) 



U{r) = Uo if Ri ^ r ^ R2, 



Uoc other- 



wise. The chemical potential ^iq — E — Uq m the disk, 
or /loo — E ~ Uoo in the leads. The rotational invari- 
ance of the problem allows us to perform the mode- 
matching for each eigenstate of the total angular mo- 
mentum = + huz/'i- (with Iz = —ihdip the orbital 
angular momentum) separately. The eigenstate of the 
Hamiltonian (21 1 corresponding to the j-th eigenvalue of 



Jz can be written as 



FIG. 4: Large graphene flake probed by two circular leads of 
radius r separated by the distance I (a) , and by a lead placed 
in the distance I' from a straight interface between undoped 
and heavily-doped regions (b). 

A similar, power-law decay of the conductance with 
the sample length is predicted for geometries with non- 
coaxial leads, considered in AppendixjX] In the two limit- 
ing situations, the Mobius transformation (All maps an 



infinite plane (hemiplane) with two (one) narrow circu- 
lar openings onto the Corbino disk. Physically, these two 
situations correspond to the setup consisting of two cir- 
cular leads probing a large graphene sample (see Fig. |4^) , 
and of one circular lead and a long straight interface be- 
tween the undoped and the heavily-doped region playing 
the role of a second lead (see Fig. BId). In the first case. 



- p»(j-i/2)v 



AJ~-J^)^^^^r), (22) 



where j is a half-odd integer j = ±i,±|,.... For the 
electron doping {E — U{r) > 0), the radial components 
X = [xAiXb]'^ for the incoming and outgoing waves are 
given (up to the normalization) by 



^.!-\/2(P) 



f 1 21 

where Hh (p) is the Hankel function of the 
(first, second) kind, and the dimensionless radial coordi- 
nate is p= \E — U{r)\r/{hvF) (so p = kr in the disk and 



xT = 



^Hfl„2{p) 



(23) 
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with 
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FIG. 5: Chemical potential dependence of the conductance 
(a) and Fano factor (b,c) at a fixed radii ratio R1/R2. Solid 
and dashed lines on panels (a,b) correspond to R1/R2 — 0.1 
and 0.5, respectively. The dotted line on panel (a) is the semi- 
classical approximation for the conductance. Panel (c) shows 
the Fano factor as a function of the chemical potential in the 
units of hvp /{R2 — Ri), with R1/R2 specified for each curve 
on the plot. 



p = Kr in the leads, with K = |/ioo The ra- 



dial current density is jr = evp ( ijj 



'i(out) 



in(out) 



Df = Im 



±H^%/,{kR^)H^'2i/2ikR2 



(25) 



Eqs. ([8|9| for G and F remain unchanged, since we 



again observe the symmetry T_j = Tj. In addition, the 
particle-hole symmetry rj(— /lo) = 2],(/io) allows us to 
limit the discussion to /uo > 0. 

Numerical values for the conductance and Fano factor 
of the doped disk are presented in Fig. [5j Following the 
idea of Kirczenow [SO. , we compare (in Fig. ^k) the exact 
quantum conductance given by Eqs. ([s]) and ( |24[ ) with the 
semiclassical approximation for large angular momenta 



Gs-cl — 



8e^ 
h 



Ji 



(26) 



±Aevp/{TTp), where the upper (lower) sign is for 
(V'™'), = ax cos If + (jySunp, and we use the iden- 
tity lin[Hl^\p)H^^_l^{p)] = 2/(7rp). For the hole dop- 

ing {E - U{r) < 0), the wavefunctions are x'''^""*^ = ^^PP^^S technique, is reproduced 



where ji = int(fci?i — 5) + 5 is the maximal value of 
j such that j ^ kRi. Surprisingly, the quantization 
steps of Gs-ci (dotted black line) are missing in the ac- 
tual data even for an extremely small radii ratio (solid 
blue and dashed red line for R1/R2 = 0.1 and 0.5, re- 
spectively). Instead, weak modulation with a period 
w TihvF / {R2 — Ri) is observed when varying /ip- Ear- 
lier, conductance quantization (with the steps of 4e^//i) 
was predicted to appear for a graphene strip with a mod- 
erate aspect ratio W / C <\ [5, .8J . The quantization with 
the steps of 8e^//i was found theoretically for a bipolar 
junction in graphene, which shows the Goos-Hanchen ef- 
fect [31]. The lack of conductance quantization observed 
here for the Corbino disk shows the role of evanescent 
modes showing a slow (power-law) decay with distance 
is crucial also far away from the Dirac point, what ex- 
hibits a striking consequence of the angular-momentum 
conservation. 

Similar to the strip geometry [7], the conductance min- 
imum at /iQ = corresponds to the maximum of the Fano 
factor (see Figs. pp). The peak width shrinks approxi- 
mately as 'Khvp/{R2 — -Ri) (for more datasets, plotted as 
a function of ^j,o{R2 — Rij/hvp = fc(i?2 — ^1), see Fig. 
|5|:). From an analytical treatment of the limit /xq — > 
for angular-momentum eigenstates (see AppendixjBj, we 
find that Eq. (16 1 for Tj, obtained within the conformal 



i(out) 



("2) 

where we use the relation Hy — 

The transmission and reflection amplitudes are obtained 
by wavefunction matching aX r — Ri and r — R2. (Note 
that the dependence of the spinor ( 22 ) plays no role 
for the mode-matching analysis.) 

Details of the calculations are given in Appendix|B] For 
|/Ltoo| 00 (the heavily-doped leads limit), the transmis- 
sion probability for the j-th mode Tj = Tj{po) reads 



C. Comparison with the Schrodinger system 

As the Corbino disk containing Dirac fermions, de- 
scribed by the Hamiltonian (21 ), has not been studied in 



16 



1 



TT^k^RlR2 (S) + )2 + (S)T)2^ 



(24) 



the literature yet, a comparison with the corresponding 
Schrodinger system is desirable for the sake of complete- 
ness. The existing theoretical works [30], however, focus 
on the model with a special, angular-momentum depen- 
dent effective potential which simplifies the analysis, but 
makes a relation to the Dirac system studied here un- 
clear. For this reason, we now present a mode-matching 
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271 4n 2n 4% 

^/2^^{R2-Ri)/h 



FIG. 6: Chemical potential dependence of the conductance 
(a) and the Fano factor (b) for the Corbino disk in a 2DEG. 
Different lines in each panel correspond to different values of 
the radii ratio: R1/R2 = 0.1 (solid blue lines), 0.2 (dashed 
red lines), and 0.1 (black dash-dotted lines). The electrostatic 
potential step is fixed at ^2m*((7o — Uao)Ri /h = 7. 



analysis for two-dimensional nonrclativistic electron gas 
(2DEG) arranged in a Corbino setup with the same po- 
tential profile as applied to Dirac fermions earlier in this 
paper. 

The Schrodinger equation for the Corbino disk in 
a 2DEG reads 

-V2 + C/(r) 



2m^ 



(27) 



where is the effective mass, and the electrostatic po- 
tential U{r) is chosen identically as in the Hamiltonian 
(21 1. The solutions are written in the form of orbital- 



momentum Iz eigenstates 



(28) 



with I integer, and the radial wavefunction <i>;(r) a com- 
plex scalar. The propagating modes in the leads exist 
only for ^loo > 0, and have the form $J"(r) = iJ^ '{Kr) 

and $°"'(r) = H[^\Kr), where K = ^2mJ^oo/^^ 
and we assume scattering from the outer lead. For the 
disk area, two linearly independent solutions are given 
by Hl^\kr) and Hf\kr) (with k = ^2m^|^o|/^^) for 
/xq > 0. Otherwise, the solutions are given by modified 
Bessel functions Ii{kr) and Ki{kr). The mode-matching 
analysis is carried out separately for each value of I [32] , 
leading to the transmission propabilities 

Ti = ( „ 1 , (29) 



where 



Ml = tI^''\k, k, R2)Tl^'^'{K, fc, i?i) 



7(2,2)/ 



T^''^\k, k,R2)Tl^'''{K, fc, i?i), (30) 



(2,1), 



^p)(X,fc,r)=[i/«(ifr) 



<i(^0 



Cl'\kr) 



- {k/K)HP{Kr) fcH\(fcr) TC}fi(fcr-) 



(31) 



with i,j = 1,2, and the upper (lower) sign corresponding 
to fiQ > {fiQ < 0). We further define 



= BMH^''^ + ^(-)'e(-Mo)ifi, (32) 

with the step function Q{x) = 1 for x > or Q{x) — 
otherwise. 

Numerical values of the conductance and the Fano fac- 
tor following from Eq. (29 ) are presented in Fig. [6] [33] for 
a large but finite value of the doping in the leads, adjusted 
such that ^J^mjJJo^^lJ^ Ri / h = 7. Both G and F are 
plotted as functions oik{R2 — Ri) ~ ^/2■m^,fJ.o{R2 — Rl)/fi 
for fixed values of the radii ratio R1/R2 = 0.1, 0.2 and 
0.5 (sohd, dashed, and dash-dot lines, respectively); G 
is additionally rescaled by a factor (i?2 — Ri)/Ri to il- 
lustrate its asymptotic behavior for kRi ^ 1, which is 
insensitive to the ratio R1/R2 [31]. We also limit the 
discussion to no > 0, as the propabilities T; given by Eq. 
(29 1 decay rapidly for /xq < 0, due to lack of propagating 



modes in the sample area. 

Although the values of G shown in Fig. [6j^a) are close 
to the semiclassical result ^34j , the quantization steps are 
absent in the data. Instead, we observe Fabry-Perot os- 
cillations with the amplitude increasing with R2/R1, for 
either G or F. We attribute the conductance quanti- 
zation, reported in earlier works |30j to the particular 
choice of the effective radial potential (note that the ex- 
isting experiments for the Corbino disk in a 2DEG [3S] 
found no conductance quantization). The main difference 
in transport through the Corbino disk geometry between 
massless fermions in graphene and massive fermions in 
a 2DEG, is the reduced backscattering at the contacts 
and the absence of details of the leads (i.e. the doping) 
in the former case [35]. This is a direct consequence 
of the energy-independent velocity in graphene, which is 
also responsible for the Klein-tunneling phenomena [371 . 
Moreover, we note a suppression of the Fabry-Perot os- 
cillations for the relativistic system. 



IV. QUANTUM BILLIARDS BOUNDED WITH 
SMOOTH EDGES 

A. Section of the disk and circular quantum dot 



A simple generalization of the formula ( 14 ) leads to the 



conformal transformation that changes a finite section 
of the Corbino disk with the inner radius the outer 
radius R2, and the spanning angle 1? (shown in Fig. [T|d) 
into a rectangle of the width W and the length £, which 
is given by 



w 

— Log — , 
§ ^Ri 



(33) 



For z = x + iy, where ^ x ^ C and < j/ ^ W, we get 
Ri ^ \w\ < i?2 and argw < ?? (with < < 27r), 
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under the condition that R2/R1 = e^'^l^ . Using 
such a condition leads to the functional K{z{w)} in the 
form 



A = A(i?i,i?2,??) 



Ri 



t/iJ 



(34) 



Thus, substituting (34 1 into ([s]) leads to the exact expres- 
sion for the conductance of a section of the Corbino disk. 
Notice that the transmission probabilities Tj for the full 
disk (16) are not reproduced for d = 27r, as the mass 



confinement is now present in the system. Instead, they 
are equal for -d — t:, what causes the conductance of such 
a half-disk to be equal to half of the full disk conductance 
for arbitrary R1/R2. The pseudodifFusive limit InA <C 1 
is reahzed for Ri « i?2, and the conductance (10 1 is 



G ~ Gdiff 



ln(i?2/i?i) 



(35) 



The above formula coincides with (17 1 for 1} = 27r. The 



opposite, quantum tunneling limit (InA 3> 1) is reached 
for Ri <^ i?2 where formula ([T3|) gives 



G « 47rcrn 



R2 



(36) 



In this case, the conductance decays (at fixed Ri) with 
the characteristic length L ^ R2 as G oc L~'"l'^ . The 
reciprocal decay, observed in Section [lTT| for the full disk, 
now appears at "d — ix. 

As a next example, we consider the conformal trans- 
formation, which changes the quantum dot shown in Fig. 
[ijc) into a rectangle of the width W and the length £. 
The transformation is given by the formula |25j 



W ^ ?« + r 

—Log , 

Woo w -r 



(37) 



with the condition (i?2 - i?i + rf l{R2 - R\ - = 



which leads to 



A(i?i,i?2,^9) 



Ri + R2 

i?l — R2 



(38) 



The origin of the coordinate system of Fig. |2] is now 
shifted to zq = {C + iW)/2. The poles of the trans- 
formation (marked by white dots in Fig. [T]d) are placed 
at w = ±r, with r = \J — R\. The angle i?oo — 
i^ooiRi, R2,'&), at which the dot edges intersect each 
other, is 



^00 = 27re(7? - i?o) - sgn(i? - ^0)^^, ^o) , 



with (.i^ji^o) — 2arcsin 



sin I sin ^ 
1 - cos f cos ^ 



(39) 



(40) 



and t?o = 2arccos (i?i/i?2)- Again, substituting (38 1 into 



n8[ provides one with the exact expression for the system 



conductance, which reaches the pseudodiffusive limit for 



Ri ~ i?2, where G ~ Gdiff [see Eq. (10)] and 



Gdif 



(41) 



In [(i?2 -Ri + rY/{R2 - Ri - r)2] ' 
whereas for the quantum tunneling limit i?i <^ R2, the 



formula (13) reads 



G fa 47rCTo 



r + Ri- R2 
r- R1 + R2 



27r/i9, 



(42) 



This leads to an asymptotic form G ^ {Ri/R2)^^^^ , as 
"i^oo ^ and r ^ R2 for R1/R2 — > 0, while the Fano 
factor approaches the Poissonian value w 1. 



B. Numerical results 



(a) 





FIG. 7: The half-Corbino disk (a) and quantum dot with 
circular edges (b) realized on a honeycomb lattice. Shadow 
areas mark heavily-doped graphene leads. Thick black lines 
indicate the mass confinement. 

We now test the analytical predictions reported earlier 
in this section by comparing them with the results of 
a computer simulation of electron transport in graphene. 
The discussion starts from the tight-binding model of gra- 
phene, with Hamiltonian 



(43) 



The hopping matrix element r^j = — r if the orbitals \i) 
and \j) are nearest neighbors on the honeycomb lattice 
(with T = 2.7 eV), otherwise r^- = 0. The single-particle 
potential Vj is arranged such that the chemical potential 
fij = Ep — Vj = /ioo in the leads marked by shadow areas 
in Fig. [tI whereas between the leads (white area) /ij = 0, 
except for the small regions, where we put fij — ha.b 
(with — — Mb> depending whether the atom belong to 
the Aox B sublattice) to model a mass-confinement on a 
honeycomb lattice [13]. Such regions are: the outermost 
edge atoms in the case of the half-Corbino disk (Fig.|7^), 
and the atoms placed out of the dot edge (thick lines 
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FIG. 8: Conductance (a,b) and Fano factor (c,d) for the half- 
Corbino disk (left) and circular quantum dot (right). Solid 
lines show the results obtained by numerical summation of Q 
and Q over the modes, dashed lines show the pseudodiffusive 
limits (10 1 and (|12[). Datapoints on left/right panels are ob- 



tained from a computer simulation of transport through the 
system of Fig. [7]; a) and[7|b). 



cial when discussing the electronic structure of closed 
quantum dots in graphene. Moreover, a recent numer- 
ical study shows that the mass confinement leads to a 
strong suppression of weak localization in such systems 
[H] . as observed earlier in experiment 42J. 

Below, we extend our numerical analysis to open sys- 
tems that cannot be obtained from a strip by conformal 
transformation, to illustrate the generic character of the 
quantum-tunneling transport in undoped graphene. 



ELECTRON TRANSPORT ACROSS A LONG 
NANORIBBON 




in Fig. ^jp) for the case of a quantum dot with circular 
edges. 

We have calculated the transmission matrix numer- 
ically by adapting the method developed by Ando for 
a square lattice [35^ to the honeycomb lattice. The re- 
sults of our computer simulation [39_, depicted by dat- 
apoints in Fig. [s] match theoretical predictions (sohd 
blue lines) as long as the number of modes in the narrow 
lead A^i > 20. Moreover, the formulas (35 1 and ( [41] ) for 
the pseudodiffusive conductance (dashed red lines in Fig. 
[8^,b) reproduce the full expression (|8| with 1% accuracy 
for i?i ^ 0.29i?2 in the case of the half-Corbino disk, and 
for Ri ^ 0.69i?2 in the case of the quantum dot with 
circular edges. Analogously, the pseudodiffusive value of 
the Fano factor F « 1/3 (see Fig. ^,d) matches the full 
expression ([9| with 1% accuracy for Ri ^ 0.43i?2 and 
i?i ^ 0.81i?2, respectively. 

In other words, the half-Corbino disk, attached to 
one narrow and one wide lead, represents the case in 
which electron transport demonstrates the pseudodiffu- 
sive character in a surprisingly wide range of the system's 
geometrical parameters. On the contrary, in the case of 
the circular quantum dot attached to two narrow leads, 
both the conductance and the shot noise show strong 
deviations from the pseudodiffusive predictions, as the 
transport is dominated by a single mode in a relatively 
wide range of parameters. The latter represents an ex- 
ample of a graphene system for which our predictions on 
quantum-tunneling transport (such as an approximately 
quadratic decay of the conductance with R2/R1) seem to 
be particularly feasible for an experimental verification, 
also because similar systems have already been fabricated 
[40j suggesting that the role of mass confinement is cru- 



FIG. 9: Nanoribbon attached perpendicularly to the semicir- 
cular (a) and rectangular (b) leads. Each system is charac- 
terized by the lead width W and the sample area length L. 
Shadow areas mark heavily-doped graphene leads. 



In this section, we present the results obtained from 
computer simulations of transport across a long nanorib- 
bon attached to the semicircular (Fig. |9^) and rectangu- 
lar (Fig. ^p) leads, which demonstrate a striking anal- 
ogy between these systems and the circular quantum dot 
studied in the previous section. 

Each of the systems in Fig. [9] is modeled by the tight- 
binding Hamiltonian ( 43 1 . The simulation parameters 



[43] are chosen to grasp the basic features of recently 
fabricated graphene nanoribbons [44' , which have zigzag 
edges and are insulating, as the weak staggered poten- 
tial placed at the ribbon edge opens a band gap in the 
electronic spectrum [IS] . A similar effect was observed in 
recent numerical studies of long nanoribbons with weak 
edge disorder [46] or irregular edges [47] , 



A. Results for an infinitely long ribbon 

We utilize the 4-terminal recursive Green's function 
algorithm [16 , which allows us to analyze directly the 
electron transport across an infinitely long nanoribbon in 
graphene. Namely, we attached two extra leads (one from 
the top and one from the bottom, not shown) to each of 
the systems in Fig.|9] that are undoped and thus contain 
the evanescent modes only. (Notice that the chemical 
potential for the outermost edge-atoms jiA.B 7^ [43,.) 
The results are shown in Fig. 10 
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FIG. 10: Conductance (a,b) and Fano factor (c,d) ob- 
tained numerically for the system of Fig. [oja) and Fig. [9|b) 
(open and solid symbols in all panels) compared with analyt- 
ical predictions (lines), (a) The pseudodiffusive conductance 
( |41[ | (solid red line) and (451 (dashed blue line). Solid and 
dashed lines in panels (b,c): the tunneling conductance (44 1 
and (47 1, and the corresponding values of the Fano factor 
F ~ 1 — Gh/Ae^. [The relation depicted by the solid line in 
the shot noise vs conductance diagram (d).] The pseudodif- 
fusive limit F = 1/3 is shown by the black dotted line (c,d). 



The conductance of a nanoribbon attached perpendic- 
ularly to circular leads (top panel in Fig. |10[ open sym- 
bols) approaches the asymptotic formula for the circular 
quantum dot ^ with Ri = W/2, R2 = L/2, ■& = tt 
ioi W K L (solid red line). For instance, a 2% agree- 
ment is reached at W/L = 0.85. This is a consequence 
of the fact that in the absence of propagating modes in 
a ribbon, most of the current flows via the central re- 
gion of the device, and the system of Fig. becomes 
effectively identical to the circular quantum dot in the 
pseudodiffusive limit, where the role of boundary condi- 
tions is negligible. For the opposite, quantum-tunneling 
limit W <^ L the corresponding formula (42) may be 
written as 



G w 2''7rCTo (W/L) 



2-r, 



with ?7 = AW/ttL, (44) 



what agrees surprisingly well with the actual data shown 
in Fig. 10 (see inset in the top panel; solid red line and 



open symbols, respectively). Such an agreement can be 
understood when looking at the current-density distribu- 
tion, shown in Fig. [TT] Even for an aspect ratio as small 
as W/L « 0.5, over 90% of the current does not leave the 
area of a circular quantum dot (bounded symbolically 
with dashed lines). 

For the case of a nanoribbon attached perpendicularly 
to rectangular leads (Fig. [9]3), the pseudo-diffusive con- 
ductance (for W ^ L) is given by [15] 



Gdiff 



aoW 



2L 



arctan ( 



L ) 



V L ) 



In 



(45) 



which is depicted in top panel of Fig. 10 (dashed blue 
line) and matches the numerical data (solid symbols) 
within 2% accuracy for W/L > 2. An identically good 
agreement with the numerics is observed for the asymp- 
totic form of the formula (45 1 Gdis ~ <Jo {W/L + I/tt), 



showing that the infinite ribbon attached perpendicularly 
to the leads has an extra (Tq/tt conductance in compar- 
ison with the rectangular geometry considered in Refs. 

[ziiHiiiiinj. 

A brief comparison between the formula (|45| and the 



generic form of the pseudodiffusive conductance (10 1 al- 
lows us to consider the functional A = A{W/L) in an 
approximate form given by 



InA 



In 



arctan 



(46) 



Subsequently, an approximate form of the quantum- 



tunneling conductance ( 13 1 for W L is 



G « 47rcrnA 



47rCToe"2 {w/Lf 



(47) 



Again, the formula (47 1 shows a surprisingly good, ap- 
proximately 10% agreement with the numerical data pre- 
sented in Fig. 10 (see the inset in the top panel; dashed 
blue line and solid symbols, respectively), suggesting 
that the power-law (approximately quadratic) decay of 
G for large i is a generic feature for transport across 
the nanoribbon, unrelated to the particular shape of the 
leads gn]. 

The numerical results for the shot-noise power are pre- 
sented in the bottom panel of Fig. [TO] The approxima- 
tive formulas ( [44| and (47 1 are substituted in the rela- 
tion F K, \ — Gh/{4e'^), which produces the analytical 
predictions depicted by solid red and dashed blue lines, 
respectively. In both cases, the agreement with numerical 
results is better than 5% when the Fano factor F > 2/3. 
An additional insight into the nature of the crossover 
from the Poissonian to the pseudodiffusive regime is pro- 
vided by an F versus G plot (see the inset). In partic- 
ular, values of F are very close to 1 — Gh/{4:e'^) even 
for a relatively large conductance G « 2e^ /h, which in- 
dicates that electron transport is governed by a single, 
valley-degenerated mode in a wide range of the geomet- 
rical parameters {W/L < 0.5 for the circular leads, and 
W/L < 1.5 for the rectangular leads). 
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FIG. 11: Current distribution in a long nanoribbon attached to circular (left) and rectangular (right) leads, as shown in Fig. 
^Si,h). Each arrow represents the average current density for a rectangle consisting of 17 x 17 unit cells. The aspect ratio of 
both systems is W/L « 0.5. Dashed lines mark symbolically the edges of the corresponding quantum dot of Fig. [7|b). 




FIG. 12: Mode-participation ratio for transport through the 
system of Fig. [o] with W/L ~ 0.3 and finite width of the un- 
doped region Woo- Three curves on each panel correspond 
to Wrx/a = 3k (squares), 3fc + 1 (circles), and 3k + 2 (trian- 
gles), with fc-integer. Top: circular leads, bottom: rectangu- 
lar leads. Lines are drawn as a guide for the eye only. 



B. Influence of armchair edges in a finite ribbon 

So far, we have analyzed the transport across an in- 
finitely long, zigzag nanoribbon attached perpendicularly 
to the leads. To find out how the results change for the 
realistic case of a long but finite nanoribbon, we consider 



now the system of Fig. |9] with a central (undoped) re- 
gion of finite width Woo ■ The system is terminated from 
the top and the bottom by armchair boundaries which 
mix valley degrees of freedom [7], so the fourfold (spin 
and valley) degeneracy of transmission eigenvalues Tj is 
expected to be replaced by the twofold (spin only) de- 
generacy. To trace the effect of armchair boundaries in 
a quantitative manner, we define the mode-participation 
ratio 



p 

J rr. 



G 



(48) 



where we assume spin-only degeneracy in summations. 
In particular, for the quantum-tunneling limit W ^ L, 
the mode-participation ratio is P„i ~ 2 if the lowest 
mode, that governs the electronic transport, has an ap- 
proximate valley degeneracy. Otherwise, in this limit 

The numerical values of the mode-participation ratio 
(48 1 are presented in Fig. 12 We took W = 80a (provid- 
ing 30 propagating modes for fi^o — t/2), L — 150\/3a 
(so W/L « 0.3) and vary Woo- The remaining parame- 
ters are identical as in the case of an infinite ribbon, stud- 
ied before. The datapoints in Fig. [T2] illustrate a smooth 
crossover from the transport dominated by a single mode 
with spin-only degeneracy {Woo ^ W L), to the situa- 
tion with full fourfold degeneracy {W L <^ Woo)- The 
details of the evolution depends on whether the width 
Woo corresponds to the metallic (Woo /a = 3fc -f 1) or to 
one of the two insulating armchair boundary conditions 
(Woo /a, = 3k, 3k + 2). In all cases, the valley degen- 
eracy is approximately restored (P^ ~ 2, marked with 
thin black line) for Wca 2L, when the role of armchair 
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edges becomes negligible, as the current is flowing pre- 
dominantly via the central area of the system (see the 
current distribution shown in Fig. 11 right panel). 



C. Implications for the experiment 





FIG. 13: These two graphene billiards both have the same 
conductance and shot noise in the pseudodifFusive regime 
W < 2L. 

For the sake of completeness, we analyze now the trans- 
port through a graphene billiard attached to two different 
leads, one narrow and semicircular, and the other wide 
and rectangular, as shown in Fig. |13[ As before, mass 
confinement (thick black lines) is applied for the edges 
not connected to the leads (shadow areas). The system 
shown in the left panel of Fig. [13] can be exactly mapped 
onto a strip (see Fig. |2| by the conformal transformation 
([37]) with the conditio n jl' + 2L - W)/{1' - 2L + W) = 
^^/^, where I' = ^4^2 _ W^. This imphcs the func- 
tional 



K{W/L) 



U ~2L+W, 



(49) 



where iDoo = 27r — 2 arcsin[/'/ (2L)] is the angle with which 
edges intersect each other in the poles of conformal trans- 
formation (see Fig.[l|;, with d = n). The pseudodiffusive 
( 10 1 and the quantum-tunneling conductance now take 



on the forms 



^^'^^ ~ In [{U + 2L- W)/{U ~2L+ W)] ' ^^^^ 



with i?oo 27r for W/L 1, and 



G « Anan 



r -2L + W 
V + 2L-W 



TTCTo 



W_ 



(51) 



where the last asymptotic expression refers to the limit 
W/L 0, for which the system behaves effectively like 
the half-Corbino disk with the inner radius Ri — W/2 
and the outer radius R2 — 2L. As a consequence, the 



pseudodiffusive conductance (50 1 reproduces the values 
obtained from the exact expression |50j with 1% accuracy 
for W/L > 0.82. 



Another striking feature of the pseudodiffusive regime 
W < 2L, which coincides with findings presented ear- 
lier in this section, is related to the fact that poles of 
the conformal transformation ( 37 1 approach the circular 
lead tip for W/L — > 2. This is why the current is flowing 
mainly through the central area of the system, and the 
two billiards shown in Fig. [13] become equivalent in such 
a limit. The earlier flndings for the circular quantum 
dot and the long nanoribbon attached perpendicularly 
to the leads, allow us to expect that the conductance of 
the nanoribbon-like system shown in the right panel of 
Fig. [13] will not deviate signiflcantly from the expression 
(SD) also in the tunneling limit W ^ L. The numerical 
results presented in Fig. [T4|conflrm such an expectation, 
as the conductance obtained by a computer simulation 
for the nanoribbon-like system |51j match again the an- 
alytical predictions for the system with circular edges 
in surprisingly wide range of the parameters. Namely, 
10% agreement is reached for W/L ^ 0.4, whereas for 
W/L > 0.8 the deviation drops below 2%. 
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FIG. 14 

shown in Fig 
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Conductance and shot noise of the two systems 
The exact (solid line) and 
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Main panel 

the pseudodiffusive (dashed line) value of conductance for a 
billiard with circular edges, and the results of computer simu- 
lation for a nanoribbon-like billiard (datapoints) . Inset: Shot 
noise vs conductance diagram for a nanoribbon-like billiard. 
Solid and dashed lines mark the asymptotic values for the 
tunneling and the pseudodiffusive limit, respectively. 



We predict that the pseudodiffusive conductance (50 1 



remains unchanged for a wide class of irregular graphene 
billiards of shapes fltting between the two limiting cases 
shown in Fig. 
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Moreover, an approximate agreement 
should be observed even for the tunneling conductance 
(51 1. We believe that such an extra flexibility in device 



setups will facilitate experiments with better agreement 
with the theory as achieved so far for rectangular samples 
|9l[T0]. In particular, the setup consisting of the narrow 
semicircular lead on one side and the straight graphene- 
lead interface on the other side, eliminates the difficulty 
of manufacturing the two parallel interfaces — one of the 
main problems that have limited the number of experi- 
mental samples, suitable for both ballistic conductance 
and shot-noise measurements, to just a few so far. 
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VI. CONCLUSIONS 



In conclusion, we have identified a novel type of quan- 
tum tunneling efi^ect, which appears in transport through 
the Corbino disk and quantum billiards in undoped 
graphene, provided that at least one of the leads (or 
billiard openings) is much narrower than the distance 
between openings L, which defines the length-scale of 
the sample. In such a tunneling limit, the conduc- 
tance G shows a slow power-law decay with L charac- 
terized by a geometry-dependent exponent. The Fano 
factor F exhibits a crossover from the pseudodiffusive 
{F — 1/3) to Poissonian {F = 1) shot noise, with a re- 
lation G ~ {1 — F)xse^ /h in a surprisingly wide range 
of L. This is because electron transport in the tunneling 
limit is effectively governed by a single mode, having the 
full spin, valley, and symplectic degeneracy (s = 8) in 
the absence of boundaries (Corbino geometry) , spin and 
valley degeneracy (s — 4) if the boundary conditions do 
not scatter valleys, or the spin-only degeneracy (s = 2) 
otherwise. In particular, for the case of a ribbon which 
contains either infinite-mass or armchair boundaries, the 
valley degeneracy is restored when armchair endings are 
shifted away from the area where the main current fiows. 
We would like to stress that the relation between G and 
F allows an experimental verification of the degeneracy 
s without referring to any geometrical parameters. 



We have explored the idea of Katsnelson and Guinea 
[15j . that transmission eigenvalues could be obtained an- 
alytically for any undoped graphene fiake of a geometry 
linked via conformal transformation to a strip, for which 
the solution is known due to Tworzydlo et al. [7]. In the 
pseudodiffusive limit, we show that the eigenvalue dis- 
tribution is affected by an arbitrary conformal transfor- 
mation only via a multiplicative prefactor, and the value 
F = 1/3 is unchanged for any closed setup, provided 
that G ^ e^/h. We test the approach for the Corbino 
disk, by comparing transmission probabilities obtained 
by a conformal mapping, and within the mode-matching 
analysis for angular momentum eigenstates. To analyze a 
crossover from the tunneling to the pseudodiffusive limit 
in a confined system, we focus on two particular billiards 
(a section of the Corbino disk and the quantum dot with 
circular edges) confined by a mass. The results of our 
numerical simulation of transport through a lattice con- 
sisting of ca. 10^ carbon atoms match the expressions for 
G and F obtained by conformal mapping. Moreover, we 
generalize the approach to obtain an approximate for- 
mula for G (which reproduces either the pseudodiffusive 
values or the tunneling-limit exponent resulting from the 
simulation) in the case of an infinite ribbon attached per- 
pendicularly to the leads — an open billiard, not linked to 
the strip via conformal transformation. 
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APPENDIX A: CONFORMAL MAPPING FOR A 
GENERIC SETUP WITH CIRCULAR 
CONTACTS 




FIG. 15: Generic circular contacts (thick lines) of radii ri and 
r2 > ri misplaced by a distance d+ > ri + r2 or d_ < r2 — ri 
(dashed line depicts the inner contact interface for the latter 
case). White dots mark poles of the transformation ( |A1[ ), thin 
lines are perpendicular to each of the interfaces (and mapped 
onto the radiant lines by the transformation). 

We consider here a generic setup, containing two cir- 
cular, but not coaxial interfaces splitting undoped and 
heavily doped graphene interfaces, as depicted schemat- 
ically in Fig. |15[ In the first case, an infinite graphene 
plane is probed by the leads of radii ri and r2 (thick solid 
lines), misplaced by the distance d+ > ri + r2 (we fur- 
ther suppose ri < and heavily doped. In the second 
case, a disk-like sample area is limited by the inner lead 
of radius ri (dashed circle) and the outer lead of radius 
r2, misplaced by d- < r2 — ri (and for d^ Q the per- 
fect Corbino geometry is restored). In both situations, 
conformal mapping onto the Corbino disk is provided by 
the Mobius transformation 



where z belongs to the disk area (with the edges radii Ri 
and i?2, ^1 ^ 1^1 ^ R2) and w belongs to the sample 
area of Fig. [TSj The real parameter (3 is adjusted such 
that 

d± = ^jrl + P^ ± ^rl + (3\ (A2) 
which leads to the useful relation 

d+d- =rl-rl. (A3) 
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(Notice that we are using one form of z{w) to describe 
the two distinct situations, in each of which only one of 
the displacements {d+,d^} is a physical parameter.) 
The explicit form of the functional A = {R2/ RiY^^ 



follows from the condition that the transformation (Al) 



always maps the first contact (of radius ri) onto the inner 
edge of the disk, whereas the second contact is mapped 
onto the outer edge. For the two situations studied here 



A(ri,r2,d±) 



xi + yi X2 ± y2 



1/2 



(A4) 



where we define the variables Xn 



Va — \J (rf+ + SadJ)'^ — (2ra)^ (with Sq, = 2a — 3, a ~ 
1, 2). The remaining one of the parameters d-}, not 
listed explicitly as an argument of A, is determined by 
Eq. ( |A3t . 

In particular, for the case of two identical circular leads 
probing a large graphene plane (see Fig. |4^) ri ~ r2 = r, 
d- = 0, and the functional 



A(r,d+) 



2r 



4^2 



2r 



4^2 



(A5) 



limit. 



00 such that I' = 
r2 = const, we find from Eq. (A4|_that 



where the approximation refers to the r <^ d 
Defining I = d+ we obtain Eq. ( 19 1 of the main text 
Similarly, taking the limit r2,d 
d+ - ri 

A « (2Z7r)i/2 for r ee n < I', what leads to E^r"|20| for 
the conductance. 



APPENDIX B: MODE-MATCHING FOR THE 
CORBINO DISK 



The radial component Xj (^) of the eigenstate ipj ( 22 1 
corresponding to the total angular momentum hj {j half- 
odd integer) and energy E can be divided into three re- 
gions. For r > i?2 (the outer lead), Xj = x}; with 



iH^^ [Kr) 



, (Bl) 



where K = \E — Uao\/fiVF with Uao — > —00 [S^, and 
the reflection coefficient. Next, for Ri < r < R2 (the 
disk area), Xj = X^' with 



II _ / H^%ikr) 



Hf\,2{kr) 



(B2) 

where k = \E — Uol/hvp. Finally, for r < Ri (the inner 
lead), Xj = xJ", "^i*^ 



(B3) 



Solving the matching conditions Xj{R2) — X^i-^^) and 
X"(^i) = xf\Ri) we find 



nKR2 



-2 (^2)^: 



H 



f\/2(Pi)+*^^^ -+1/2(^1) ' (B4) 



(1) 



Here, we derive the transmission and reflection ampli- 
tudes for scattering eigenstates of the Hamiltonian (21) 
for the Corbino setup, as shown in Fig. [ija). Without 
loss of generality, we suppose electron doping in the leads 
E > Uoo, but an arbitrary doping in the sample area 
ri = sgn{E - C/o) = ±1. 



i?l\/2(Pl)+*^^+l/2(Pl) 



r(i) 



j+1/2 
(1) 



^iil/2('^l)+*^^)+l/2(/^l) 



(B5) 



where pa = kRa (with a = 1,2), K2 ~ KR2 — T^j/'^-, and 
we have deflned 



Im 



HY\,2{pi)Hf^„2{P2) ± H^,/2^p,)H^,f^p2) 
= -Jj-l/2{Pl)Yj^i/2{P2) + Yj_i/2ipi)Jj^l/2{P2) T Jj + l/2{Pl)Yj±i/2{P2) ± ^j + 1/2 (Pl) '/j±l/2 (^2) , (B6) 



(1) 



r(2) 



with Jiy{p) [Yi,{p)] the Besscl functions of the first [second] kind. The refiection and transmission amplitudes are 



X {Hj!.\/2(P2) [H^%iPi) +^vH^Z/2iPl)] - H^j%iP2) [H^%iPi)+^vH^%2iPi) 



7/1)7 - 



(B7) 



r 



and 



Defining Tj = |tj(£')r, we obtain Eq. (24) of the mam 



Arje 



iK{Ri-B.2) 



(B8) 
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text. Notice that Tj depends solely on = E — Uq, 
as tj{E) is affected by Hoo = E — Uoo only via a phase 
factor. It is also insensitive to the doping sign 77 = ±1, 
what corresponds to the particle-hole symmetry. 

For the undoped-disk limit (fc ^ 0), Eq. (B6I leads to 
the asymptotic form 



Substituting the above into Eq. (B8) we obtain Tj = 
\tj{E — > J7o)P as given by Eq. (I6I of the main text. 
Hence, the correspondence between the mode match- 
ing for angular-momentum eigenstates and the conformal 
mapping technique for the disk in undoped graphene is 
established. 



2^ 

Trk\/RiR2 



i?2 



R2 

Ri 



(B9) 
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